Chirality and Z2 vortices in a Heisenberg spin model on the kagome lattice 
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The phase diagram of the classical Ji — J2 model on the kagome lattice is investigated using exten- 
sive Monte Carlo simulations. In a realistic range of parameters, this model has a low-temperature 
chiral-ordered phase without long-range spin order. We show that the critical transition marking 
the destruction of chiral order is preempted by the first order proliferation of Z2 point defects. The 
core energy of these vortices appears to vanish when approaching the T — phase boundary, where 
both Z2 defects and gapless magnons contribute to disordering the system at very low temperature. 
This situation might be typical of a large class of frustrated magnets. Possible relevance for real 
materials is also discussed. 

PACS numbers: 75.10.-b, 75.10.Hk, 75.30.-m, 75.40.-s, 75.40.Mg, 75.50.Ee, 75.90.-(-w 



In classical spin systems, competing interactions com- 
monly frustrate the conventional (tt, tt) Neel order, pos- 
sibly leading to more exotic arrangements of the lo- 
cal spins. Prominent examples include helicoidal con- 
figurations, [TJ [21 13] which usually break space inver- 
sion and time reversal. Noticeably, such spin chirality 
is a sufficient condition for multiferroic behavior, i.e. 
non-zero coupling between magnetization and electric 
polarization [H [5]. 

However it is not uncommon that non-planar magnetic 
orders relieve the frustration even more effectively than 
helicoidal configurations do. The associated magnetic or- 
der parameter is then three-dimensional, hence also chi- 
ral. To date, two such orders have been exhibited, both 
on triangular-based lattices with competing interactions: 
a 4-sublattice tetrahedral order was found on the trian- 
gular lattice ID while a 12-sublattice cuboctahedral order 
was found more recently on the kagome lattice. [7] 

In two dimensions, complex magnetic orders might 
seem of purely theoretical interest since gapless spin- 
waves destroy the spin long-range order at arbitrarily low 
temperatures. However, this disordering process is soft, 
in the sense that at low but finite temperature, the spin- 
spin correlations remain large enough to sustain emergent 
long range orders. This is examplified by the two above- 
mentioned models where chiral long-range order persists 
up to finite temperatures, whereas long-range order in 
the spins is lost. Interestingly, the emergent chiral order 
parameter is Ising-like, and in a straightforward extrapo- 
lation one expects that these chiral phases will disappear 
through a critical transition in the 2D Ising universality 
class. However, Momoi and Kubo showed that in the case 
of the tetrahedral order, this is true only in the "weak 
universality" sense. [HI 

We point out that such three-dimensional magnetic or- 
ders completely break the SO{3) symmetry of Heisenberg 
interactions. Hence the order parameter space is 5*0(3) 
which supports point defects, namely vortices in two di- 



mensions, around which the order parameter is rotated 
by 27r. However, note that rotation of 47r is equivalent to 
the identity, so that the order parameter may only wind 
by ±27r, as can be more formally deduced from Hi = Z2. 
This evidences the pecuHar topology of 50(3) vortices 
compared to the well known 50(2) vortices involved in 
the Berezinski.-Kosterlitz-Thouless (BKT) transition In 
particular, since 50(3) rotations of 47r are equivalent 
to the identity, 50(3) only supports vortices with unit 
"winding number" , as can be deduced more formally 
through Hi (50(3)) = Z2. These Z2 vortices were first 
exhibited in an early numerical work by Kawamura and 
Miyashita [5] on the antiferromagnetic Heisenberg model 
on the triangular lattice, where the defects were shown 
to proliferate rather abruptly at finite temperature. To 
date, however, there is still no conclusive evidence that a 
genuine phase transition indeed takes place in this model. 
On the experimental front, the proof of existence of the 
Z2 vortices remains rather elusive, although they may 
have been probed indirectly in recent NMR experiments 
on NaCr203. ^ 

In this letter we exhibit Z2 vortices in the Ji — J2 model 
on the kagome lattice and show that they are responsible 
for the first order nature of the chiral transition and we 
study the effects of frustration on the core energy of these 
defects. 

The Hamiltonian of this model reads: 
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where the first sum runs over pairs of nearest neighbors 
(at distance 1 on the kagome lattice) and the second sum 
over pairs of second nearest neighbors (at distance %/3) • 
We are interested in the 12-sublattice antifer- 
romagnetic ground state obtained for Ji < and 
J2 > I'^il/S. [7> The associated order parameter has 
the symmetry of a cuboctahedron, with scalar chirality 
Cijk = V^ Si ■ Sj X 5fc either -1-1 or —1, where {i,j,k) 
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FIG. 1: (Color online) 12-sublattice Neel order on the kagome 
lattice. Top: The 12-site magnetic unit cell: different sublat- 
tices are indicated by different numbers. Bottom: the order 
parameter in spin space and its mirror image. The two degen- 
erate ground states differ only by their scalar chirality, namely 
the triple product of 3 spins on a (shaded) triangle. 



label the 3 sites of a triangle clockwise (Fig. IT]). On the 
kagome lattice, triangles pointing up and those pointing 
down carry opposite chirality in the ground state and the 
associated chiral long-range order is evidenced by a finite 
value of the staggered chirality 
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where the sum runs over all 2 A^/3 triangles of the kagome 
lattice with aijk = 0(1) for triangles pointing up (down). 
To investigate the finite-temperature phase diagram of 
this model, we perform Monte Carlo simulations using 
a parallel tempering algorithm. This method is indeed 
suitable to overcome the free energy barriers encountered 
at first-order phase transitions, as we will show is the 
case for the chiral transition in this model. Further, 
this algorithm is easy to parallelize and yields thermo- 
dynamic quantities over a large range of temperatures 
in a single run, once combined with re weighting meth- 
ods. We simulated samples of linear size L ranging from 
12 to 64, with up to TV = 3L'^ < 12288 spins. Al- 
though the tempering method suppresses the slowing- 
down associated to the crossing of free energy barriers, 
the existence of a large spin-spin correlation length drives 
an "effective critical" slowing-down. Therefore, for the 
largest samples the number of Monte Carlo steps needs 
to be increased up to 2^^ steps per spin. The first re- 
sults of the simulations are collected in Fig. |2] Start- 
ing from J2/I Ji| = 1/3, the temperature extent of the 
chiral phase increases with increasing antiferromagnetic 
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FIG. 2: (Color online) Transition temperature T* fronr chi- 
ral order to full disorder versus the coupling ratio J2/IJ1I . 
For < J2/IJ1I < 1/3, the ground-state is ferromagnetic. 
For J2/IJ1I > 1/3, the ground state is the 12-sublattice Neel 
ordered phase described in the text and abbreviated here as 
"cuboc". Dots: results of Monte Carlo simulations (size ef- 
fects are much smaller than the size of the symbols). The 
solid line is a guide to the eye. 



Ji. We have monitored several thermodynamic quanti- 
ties, namely the averaged value of the energy per spin 
the specific heat, C^/Zcs — ^.^^1^ ((^^) ~ (s)^)' ^he 



staggered chirality (|C|) and the associated susceptibil- 
ity fcsxc = If ((C^) - (|C|)2). Figure g shows the 
rapid destruction of the chiral long-range order at the 
transition (note the very small temperature scale) for 
J2/|Ji| = 0.38. 

To characterize the chiral transition we use standard 
finite-size scaling analysis. For 1/3 < Jil\J\\ < 0.45, 
the energy distribution becomes bimodal in the neigh- 
borhood of the transition, and its maxima become more 
pronounced with increasing sample sizes. Both the 
maximum of the specific heat C^^^{L) and of the chi- 
ral susceptibility x^^'^^{L) increase algebraically with L 
with exponents within 2.0 ± 0.15 for J2/IJ1I < 0.38. 
For J2/IJ1I = 0.39 and 0.40, the scaling regime is 
reached only for the largest samples. In addition, de- 
noting T(7„ (L) and T^^ (L) the temperatures of the max- 
imum of Cy{L) and Xc{L) one obtains l/Tc„^g(oo) — 
l/Tc^,xciL) ^ 1/i^ for J2/IJ1I < 0.38. The alDOve anal- 
ysis shows that the transition does not belong to the ex- 
pected Ising universality class, but is first-order, although 
the increasing difficulty to reach the scaling regime, when 
increasing J2/I Ji| , evidences the concomitant growth of 
a correlation length. This is consistent with the contin- 
uous decrease of the latent heat with increasing J2/I Ji| 
(0.034,0.028,0.012 for J2/IJ1I = 0.36,0.38,0.45) 

To account for the discontinuous nature of the chiral 
transition a thermally activated mechanism was looked 
for: proliferation of point defects is an obvious candi- 
date, and we now proceed to compute the Z2 vortices of 
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FIG. 3: (Color online) Top: Temperature dependence of the 
staggered chirality {\C\) and the specific heat for J2/IJ1I = 
0.38 and L — 64. Temperatures are measured relatively to T* , 
temperature where C'v is maximum. Bottom: spatially aver- 
aged spin-spin correlation length fspin and the vortex density 
versus temperature. 



the model. We define the local trihedron of a 12- sites 
magnetic unit cell as follows: ei = {Si + Sj)/\Si + Sj\, 
63 = {Si — Sj)/\Si — Sj\ and 6*2 = 63 x ei, with i and j 
any two sites in the magnetic unit cell that carry non- 
coUinear spins. These trihedra live on a triangular super- 
lattice of spacing 4 in units of the original kagome lat- 
tice. Then we determine the rotation of axis n and angle 
9 between two successive trihedra, and disambiguate it 
from {—ft, 2tt — 9) by using its universal covering element 
U{n,9) = e~*5'5'" e SU{2), where a are the three Pauh 
matrices. Finally we define the Z2 vorticity inside an 
elementary (triangular) contour by 
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where T^a = 1(0) when the loop A encloses a singu- 
larity (otherwise). We computed the vortex density 
{riy) = Q/N"^^ Va, where the sum runs over all A^/6 
triangles of the super-lattice. The results are displayed 
in Fig.|3]versus {T-T*)/T*, where T*{L) is the temper- 
ature of the maximum of Cy {L) . They clearly show that 
the chiral transition is concomitant with the proliferation 
of vortices. 

The behavior of the spin-spin correlation length ^spin 
close to the transition evidences that vortex proliferation 



drastically decreases the magnetic short range order as 
well as it kills the emergent chiral order, far before the 
expected critical regime of the chiral phase is attained. 
Hence the Ising chiral transition is avoided simply be- 
cause the vortex proliferation triggers a first-order phase 
transition that preempts the critical regime. We em- 
phasize that the disordering effect of the Z2 defects is 
much stronger than that of Z vortices at the BKT tran- 
sition. [TD] Fig. 3 indeed shows that ^spin is divided by 
three in a temperature range ~ 10""^ T*, and this brutal 
decrease may even be smoothed by the finite size of the 
sample. 

The chiral correlation length ^chiral i ^^ computed from 
the structure factor of the chirality C at the transition, 
measures the discontinuity of the chiral transition. For 
J2/\Ji\ — 0.38 we find ^chiral — 6, while it exceeds the 
largest available lattice size (^chiral > 64) already for 
J2/\Ji\ = 0.40. Consistently, the vortex proliferation 
smoothes upon increasing J2/|'^i| (Fig- ^ and for large 
enough J2/IJ1I = 0.45, the number of vortices at the 
transition is clearly seen to decrease. Although this is 
not the core of our study, note that the fast growth of 
Cchirai with increasing J2/IJ1I makes it particularly deli- 
cate to idenitfy a possible critical end point to the line of 
first order chiral transitions shown in Figure |2] 

Discussion: The chiral transition in the present Ji — J2 
model may be typical of many frustrated magnets. In- 
deed, as emphasized above, as long as the ground state 
completely breaks the 0(3) symmetry of the Heisenberg 
Hamiltonian, as is the case for non-planar Neel orders, 
both chirality and Z2 vortices exist. The complete break- 
ing of S0{3) induces Z2 point defects, while the breaking 
of time reversal (the discrete part of 0(3)) leads to chiral 
degenerate ground-states. Hence we expect very simi- 
lar features for the chiral transition in the present model 
and that observed on the triangular lattice by Momoi et 
al. : [^ in a posterior study, these authors indeed noticed 
that the chiral transition may be weakly discontinuous. 

However, the nature of the chiral transition cannot be 
deduced from symmetry arguments alone and will ulti- 
mately depend on the energetics of the two competing 
mechanisms that suppress chiral order: formation of chi- 
ral domain walls versus creation of Z2 vortices. 

In this respect, the case of the Ji — J3 model on the 
square lattice is interesting. Capriotti and Sachdev [TT] 
have shown that in this model, the doubly degenerate he- 
licoidal ground-state results in a finite temperature chiral 
phase. Although Z2 vortices are allowed by symmetry, 
the chiral transition at finite temperature clearly falls 
in the 2D-Ising universality class. [TT] Consistently, we 
understand that in a problem dominated by antiferro- 
magnetic long- wavelength fiuctuations, forming a chiral 
domain wall is much less costly than creating a Z2 vortex. 

Correspondingly, in the Ji — J2 model under study, 
we expect that competing ferromagnetic and antiferro- 
magnetic interactions will enhance the short-range fiuc- 
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FIG. 4; (Color online) Evolution of the vortex density {n-v} 
around the transition temperature T* as a function of J2 / 1 Ji | . 



tuations, hereby decreasing the core energy of Z2 vor- 
tices. This qualitative argument is supported by the ob- 
servation that the discontinuity of the chiral transition is 
maximal when the competition of interactions is highest, 
i.e. close to J2/\Ji\ = 1/3: there the chiral transition is 
clearly triggered by the proliferation of defects. 

Note also that the decrease of the transition temper- 
ature as J2/IJ1I -^ 1/3 suggests that the core energy of 
the Z2 vortices vanishes at the T = phase boundary, 
I.e. J2/IJ1I = 1/3. 

However, in this very low temperature regime quan- 
tum fluctuations are expected to play a significant role. 
The r = quantum problem was actually studied pre- 
viously by three of us: [7] using exact diagonalizations 
for spins 1/2 as well as spin-wave expansions, it was 
shown that both classical ground states, ferromagnetic 
and 12-sublattice antiferromagnetic, survive quantum 
fluctuations with a quantum phase transition located at 
J2/IJ1I ~ 1/3. Hence, although the present study is 
purely classical, it yields the promising prospect that ad- 
ditional excitations become gapless exactly at the quan- 
tum phase transition, supplementing the usual gapless 
magnons that exist in either of two neighboring phases. 
This is highly reminiscent of the breakdown of the Lan- 
dau paradigm in quantum spin systems, |12j although 
the route from classical Z2 vortices to some fractional- 
ized "spinon" excitations for spins 1/2 remains a totally 
open problem. 

Experimentally, a large number of magnets on the 
kagome lattice have been synthetized up to now, and it 
has been a long route to arrive at the Herbertsmithite 
Z„CuO, which remains an antiferromagnetic spin liquid 
down to 50 niK, |13( I14j with a dominant antiferromag- 
netic interaction of about 190 K. Amongst other difficul- 
ties, the sign of the nearest neighbor coupling is a pending 
problem: Atacamite, parent both of Paratacamite and 



Herbertsmithite, becomes ferromagnetic around 10 K. 
Similarly, Cutitmb, the organic compound at the ori- 
gin of our interest in the present model, [15j [16] was 
recently shown to experience a transition to three di- 
mensional ferromagnetic order around 500 mK |17j : this 
tendency to ferromagnetism is deeply rooted in the geom- 
etry of the exchange paths between nearest neighbor Cu 
ions on the kagome lattice. If the malediction of such 
low-temperature ferromagnetic orderings can somehow 
be avoided, possibly through a decrease of the inter-layer 
couplings, the present study would be of direct experi- 
mental interest, beyond its initial theoretical motivation. 

In this letter, we showed that in the Ji — J2 model 
on the kagome lattice, chiral order is wiped out at finite 
temperature by the first-order proliferation of Z2 vor- 
tices. In the region of extreme frustration J2/I Ji| — 1/3 
the core energy of the defects decreases and appears to 
vanish exactly at the T = phase boundary. This be- 
havior is probably common to frustrated spin systems in 
which competing interactions lead to three-dimensional 
antiferromagnetic order parameters. 
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